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A left-unilateral matrix equation is an algebraic equation of the 

form 

O-i. flo + aix + a2X^ + . . . + a„x" = 

D , 
F^ ■ where the coefficients a^ and the unknown x are square matrices of 

the same order and all coefficients are on the left (similarly for a right- 

k>( \ unilateral equation) . Recently certain perturbative solutions of unilat- 

;_i ' eral equations and their properties have been discussed. We present a 






> 
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unified approach based on the generalized Bezout theorem for matrix 
polynomials. Two equations discussed in the literature, their pertur- 
bative solutions and the relation between them are described. More 
abstractly, the coefficients and the unknown can be taken as elements 
of an associative, but possibly noncommutative, algebra. 
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1 Introduction 

In a discussion of generalized Born-Infeld theories |I|, the construction of 
the Lagrangian was reduced to the solution of certain unilateral algebraic 
matrix equations and it was conjectured that the iterative solution of those 
equations is a sum of symmetric polynomials in the coefficients and of terms 
which are commutators. Equivalently, that the trace of the matrix solution is 
equal to the sum of traces of symmetric polynomials in the coefficients. The 
conjecture was later proven in and by A. Schwarz in [Q using different 
methods and a slightly different form of the equation. 

In the present note we combine Schwarz' s idea of expressing the trace 
of the solution as a contour integral in the complex plane of the trace of the 
resolvent of the corresponding matrix with the idea used in |^ of using the 
basic property of the trace of the logarithm of matrices. 

The two approaches can be easily combined by using the generalized Be- 
zout theorem for matrix polynomials. The qualitative algebraic fact that the 
trace of the solution is given by a sum of traces of symmetrized polynomials 
in the coefficients emerges almost without any computation. The coefficients 
in the expansion are not hard to compute and we give an explicit expres- 
sion for them. Q The method used here applies equally to the two different 
equations considered in and 0]; in the last section we clarify the relation 
between the two equations and their solutions. 

Before closing this section we recall the statement of the generalized Be- 
zout theorem (see e.g. p). Let 

P(A) = ao + aiA + aaA^ + . . . + a„A", \eC (1.1) 

be a polynomial with square matrix coefficients a,, and x a square matrix of 
the same order. Define 

P{x) = Co + aix + a2x'^ + . . . + a„x" (1.2) 

with the coefficients all on the left. It is easy to verify that 

P{X)-P{x) = Q{X,x){X-x), (1.3) 

where 

n— 1 / '^ \ 

Q(A,x) = J2^\ Yl e«ra^"~'"^ • (1.4) 



/=0 \r=Z+l 



^Our contour integral formulas are very similar to Schwarz' s, but there seem to be 
some minor discrepancies between ours and his. 



In other words A — x is a divisor of P{X) — P{x) on the right (if we had taken 
all the coefficients on the right, then X — x would have been a divisor on the 
left). 

We shall study matrix equations of the type 

P{x) = , (1.5) 

which we shall call left-unilateral matrix equations (meaning that all the 
coefficients are on the left). If x is a solution of (|1.5|) , the characteristic 
polynomial P{X) of the equation can be factorized as 

P{X) = Q{X,x){X-x). (1.6) 

It is clear that the Bezout theorem applies more abstractly if one considers a^ 
and X as elements of an associative, but possibly noncommutative, algebra. 
The same remark applies to the rest of this note, if one uses an appropriate 
algebraic definition of the trace as cyclic average (see |0]). 

2 The Trace of the Solutions of Unilateral 
Matrix Equations 

A. Schwarz |^ considers the unilateral matrix equation 

x" = 1 + e (ao + aix + . . . + a„_ia;""^) . (2.1) 

where e is a small parameter. For e = the equation has n solutions in the 
complex plane, the n roots of 1. For e 7^ each of these solutions admits a 
perturbative expansion as a formal power series in e. A. Schwarz proves that 
Tr X* can be expressed as a power series whose coefficients are symmetrized 
products of the Oj and gives an expression for the coefficients in terms of 
contour integrals in the complex plane. 

We have found that a more explicit expression for Tr x^ is provided by 

°° e'^ TriSfa"" a"""^) ^~^ ( """"^ 

Tr x"* = Tr 1 + s > —r > ; ^ — TT s + > Irii — rn 

A^ n'' ^-^ no!...n„_i! J-i I A^ ' 

k=\ no+---+nn-i=fc r=\ \ l=\ 

(2.2) 
This formula holds for positive as well as for negative values of s. Here the 
expansion is made around 1. A similar expansion could be derived for all the 



n iterative solutions of the equation. The normahzation of the symmetrized 
product iS(aQ° . . . a"") is chosen in such a way as to give the ordinary product 
if the factors commute P . 

If we introduce the notation 



n^^^^^ for;^ = l,2, 



r=l 



(2.3) 



for A; = 



which reduces to the usual definition of ( , 1 for a G N, then the result (p.2|) 
can be rewritten 

Trx" = 



(2.4) 



Tri+'-Y^e\k-i)\ y: 



Tr5(a^° . . . <""/) / n(^ - '^ + Er=i ^^i 



k=l 



no!...n„_i! 



?iO+...+n„_i=A; '^ " ^ \ k — 1 

As a first step to prove ( p.2[ ), we apply the generalized Bezout theorem. 

(2.5) 
(2.6) 



For (|2.1| ) the characteristic polynomial is 

P{\) = 1 - A" + e (ao + aiA + . . . + a„_iA"'^) 
and, if Pix) = 0, 



P(A) = Q(A,a;)(A-x) 



with 



n— 1 / n— 1 



Q(A, x) = Y,\ Yl eflrx^'"'-^ - x"-^-i A' . (2.7) 

1=0 \r=l+l J 

Let X be one of the solutions of ( p.l|) , e.g. the one which reduces to 1 for 
vanishing e. (For the other solutions a similar procedure can be followed). 
We take the logarithm and then the trace of ( p.6|) , and obtain 

(2.8) 



TrlogF(A) = Trlog(Q(A,2:)) +Trlog(A-x) . 
Differentiation with respect to A yields 



(2.9) 



More generally, we can multiply the above equation by any function /(A) of 
A, which is regular in a neighbourhood of 1, so as to get 

Tr^P'(A)/(A) = T.^Q'(A..)/(A) + T.^/(A) . (2.10) 

As a next step, we follow Schwarz' s idea of using a contour integration 
to isolate the relevant part in the trace. We consider a small circle F around 
1 in the complex plane, or more generally a small closed curve winding once 
around 1, and containing no other solution of ( p.l|) for e = 0. Then we 
compute the integral of equation ( p.lO|) along it for small e 



i2my' j( dX Tr-^P'(A)/(A) = (2.11) 

{2m)-' j( dX Tr— i-yQ'(A, x)/(A) + (2^)-' £ dX Tr^/(A) . 

The two integrals on the right hand-side can be evaluated through the Cauchy 
theorem in the following way. We are considering the case of small e, and then 
Q is a polynomial in x and A, which vanishes for A near e~^ , k = 1, . . . ,n — l, 
so that in this case its inverse Q~' has no singularities near 1. Q' and / are 
regular functions near 1 and therefore the Cauchy theorem guarantees that 
the first term on the right hand-side vanishes: 

{2m)-' I dX Tr— ^Q'(A, x)/(A) = . (2.12) 

On the other hand Tri^/(A) has poles for A near 1 with total residue Tr/(x). 
As X is close to 1 for small e, the Cauchy theorem yields 



{21X1)-' i dXTi- /(A) = Tr/(x) . (2.13) 

Jv A-x 



Finally, we obtain 



Tr/(a;) = {2m)-' j( dX Tr-^P'(A)/(A) . (2.14) 

Therefore, the problem of computing Tr/(x) amounts to evaluating the in- 
tegral on the right hand-side of ( |2.14| ). 



We can factorize 

P(A) = (1 - A")r(A) , (2.15) 

where 

n-l 

r(A) = l-e(A"-l)-^^azA'. (2.16) 

1=0 

Then 

1 nA""^ 1 

and ( [2.14] ) becomes 

The first integral on the right hand- side equals Tr /(I). To evaluate the 
second term we integrate by parts 

±.£dXTT^^r{X)f{X) = -^ j((iATrlog(T(A))/'(A) . (2.19) 

This is justified, because TrlogT(A) and /(A) are regular on F, so that 

-i, / dA Tr-^ (log(T(A))/(A)) = 0. (2.20) 

Im Jy dX 

In this way we obtain 

Tr f{x) = Tr /(I) - ^ (f dXTi log(T(A))/'(A) . (2.21) 

We expand the logarithm 

OO h. 

— for a small (2.22) 

k 



k=l 



and then make the following change of variable 

2/ = A" . (2.23) 



This is possible, because if we start with a small closed curve T winding once 
around 1, then after the variable transformation (|2.23|) we still have a closed 



curve winding once around 1. We restrict ourselves to the case /(A) = A* 

^ ZdATrlogfl-e^^kL^l A^-i= (2.24) 



2m Jy \ A" - 1 

n{^^2mkJr \ y-l J ' 

Notice that y^ is regular in a neighbourhood of 1. In this form we can already 
see that the result is symmetrized in the coefficients {at}, because they enter 
through expressions of the type 

, (=0 / no+...+n„-i=k 

We obtain 
Trx" (2.26) 



^ n—l 

Tri + ^Ve'^ y ^'^7 ---V^ :. . y '=^ 



71 



n -^^ ^^-^ uqI . . . n„_i! rfw'' 

fc=l no + ...+nn-i=k 



y=l 



oo git TriSfa"" a""~^) '^""'^ / "~"^ 
= Tr 1 + s > — r > r ^^T^ IT -5 + / ^"-« ~ '""' 

fc=l no + -..+n„_i=fc r=l \ 1=1 

Here we have applied the Cauchy theorem in its more general form 



c 



(y - Z/o)^ k-l\ dy^ 



where C a closed curve winding once around y^, and f{y) is a function which 
is regular inside C. 

6 



In the simplest case n = 2 it is possible to solve the classical equation 
(i.e. the equation for commuting aj) and this provides a closed expression 



for (O) 



Tr X = Tr 



e| + 5jl + eao + (e|)2 



(2.28) 



where we again have chosen the solution which reduces to 1 for e = 0. In 
the appendix we show how expanding the square root in (|2.28|) near 1 we 
actually recover ( |2.2| ) in this case. 

3 Relation between two unilateral matrix 
equations 



The same procedure which we have used to prove ( |2.2| ), namely applying the 
generalized Bezout theorem and then the Cauchy theorem, can be equally 
well followed for the equation 

<l> = Ao + A<l> + ...A„$", (3.1) 

which is studied in P] . If we define 

A{\) = AQ + A^\ + ...AnX' (3.2) 



the characteristic polynomial of ( p.l|) is A — A[X) and the result which cor- 
responds to ( |2.21| ) is 



Tr/(<|.) = -^.i rfATrlog(l - ^)/'(A) . (3.3) 

Here C is a closed curve winding once around 0. We consider the particular 
case /(A) = A** and expand the logarithm 

Tr $^ = , y i^:ill: y Tr5(y...yl^) f ^^ ^^Er^oC-iK+^i . 

k=\ no+...n„=fc 

(3.4) 
The Cauchy theorem has the effect of selecting the words of dimension 
^JLq(^ — l)n; = —s in the expansion of the logarithm, i.e. 

oo ^ 

Tr $^ = .Tr J] -(Ao + . . . A„)lELoa-iM=- (3-5) 

fc=i 



and we recover the result already stated in 0, where the concept of dimension 
of a word was introduced. It should be remarked that the solution to ( |3.1|) 
for vanishing coefficients is $ = 0. Consistent with that, as already noted in 
in 1^, the result (|3.5|) only holds for positive integers s: it has no sense to 
invert the solution in this case, but see below. 

To study more closely the relation between equation ( |2.1|) and (|3.1| ) we 
make the Ansatz 

x = l + a$. (3.6) 

(Again, the same procedure could be followed for any of the roots of unity 
which are the solutions of ( p.l|) for e = 0.) 

We choose the parameter a in such a way that 

a"-i = -n . (3.7) 

There is no reason for a to be real, since most of the solutions of ( p.l| ) are 
not real, even for real coefficients. 

Then the relation between the n coefficients Oj and the n + 1 coefficients 
Ai is easily found to be 



Ai 



n~l 



-a'-"e 



y [ J ] o-r for / = 0, 1 



n-1 

n \ ,_„ sr^ I r 



r=l 



a' " \ , I — a' "e\^ I . j a^ for 2 < / < n . 



(3.8) 



Some remarks can be made with respect to ( |3.8| ). First, observe that An = 1 
is fixed, but this had to be expected, because the equation ( ^.Vj ) has one 
coefficient less than the equation (pTTp. 

Moreover, (|3.8|) is a linear relation, and it is invertible, so that once the 
result that the trace of the powers of a solution depends only on the sym- 
metrized products of the coefficients is proven for one of the two equations, it 
immediately follows also for the other. Negative powers of x can be expanded 
from (|3.6|) into a series of positive powers of $. 

If e is small, only the first two coefficients Aq and Ai are automatically 
small, the other coefficients satisfy 

Ai '^° «'-" ( "}] for 2 < / < n . (3.9) 



Therefore, the two expansions ( |2.2| ) and ( |3.5| ) do not necessarily hold for the 
same range of the coefficients. 
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Appendix 

In this appendix we check that (|2.2| ) actually is the series expansion of (|2.28 



for s = 1, n = 2. We start by expanding the square root in (|2.28|) around 1 



Tr a; = Tr 



r=l 



:Trl 



Tr ai 



(A.l) 



oo _ r 



r2r-no 



no!(r-no)! 



r=l no=0 

We introduce the new variables: 

rii = 2(r — no), k = uq + rii = 2r — uq . 

Then we can rewrite 

^ ^ Trai 
Tr X = Tr 1 + e 



Tr^la^Ml)'^'"" 



(A.2) 



(A.3) 



k=l no+ni=k 

nieven 



-1^ 



^+no-lK-l)" 



no!ni! 



(2no + ni-3)!!Tr5(a;j°ar), 



where we have used the relation 

rill 



f!2^ 



(rii — 1)!! for rii even. 



(A.4) 



Due to the conditions no + rii = k, rii even, the following relation holds 

fe-i 
(_l)^+"o-i(„^ - l)!!(2no + ni- 3)!! = JJ (1 + m - 2r) (A.5) 

r=l 

and so formula ( |A.3| ) can be brought into the more compact form 

OO L. k-1 

Trx = Trl + ^- J^ ^Tr 5(a^»ar) J] (^ + ^i " 20 , (A.6) 



2'^ ^ .rioini! 

fe=l no+ni=K 



r=l 



which now coincides with (p.2| ) for s = 1, n = 2. It is no longer necessary to 
explicitly sum only over even values of ni, because 

fc-i 

JJ (1 + rzi - 2r) = for r^i odd, 1 < ni < k, k > 1 . (A.7) 



r=l 
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